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A model for the pd — > pdr] reaction published earlier, including the final state in- 
teraction (FSI) of all particles, is revisited to investigate the low energy data on this 
reaction. The three body problem of p-d-r] scattering in the final state is approximated 
in terms of pairwise interactions between the three particles in the final state. Apart 
from a comparison with some preliminary data, two new findings relevant to the near 
threshold data analysis are reported. The first one points toward the limitations of an 
FSI factor used conventionally to extract the eta-deuteron scattering length and infer 
subsequently on the existence of eta-mesic states. The second result emphasizes the role 
of the p — d FSI and the strong Coulomb repulsion near threshold. Finally, a comparison 
of the above model calculation with low energy data, excludes very large eta-nucleon 
scattering lengths. 
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1. ri meson production 

The strong and attractive nature of the rj — N interaction in the s— wave^ mani- 
fests itsel f in the sh arp enhancement of cross sections in yy-producing reactions near 
threshold l ^ l ^ | 4 | 5 | 6 | 'pj^jg enhancement was attributed to the possible existence of 
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exotic ry— nucleus quasi-bound states ' first proposed by Q. Haider and L. C. Liu. 
The experimental findings motivated theoretical searches of metastable states of 
77-mesic nuclei (see Refs. 7-15 for some theoretical predictions and an extensive list 
of references of other theoretical searches) , as well as experimental searches for the 
evidence of such states in the eta meson producing reactions. 

To investigate this possibility in the three nucleon system, in the recent past 
total as well as differential cross sections for the pd pdrj reaction have been 
measured'-" at the CELSIUS storage ring of The Svedberg Laboratory, Uppsala, 
using the WAS A/PRO MICE experimental apparatus for excess energies (Q) > 14 
MeV. Due to the strong energy dependence of the total cross section close to thresh- 
old, the data at low energies are however better suited to study the effect of the 
final state interaction (FSI) and explore the existence of quasi-bound eta-nuclear 
states. For Q < 14 MeV, two data sets exist for the pd pdrj reaction. One 
measurement was performed with the SPESIII spectrometer at SATURNE for two 
excess energies, namely, Q — 1.1 MeV and 3.3 MeV ^1^. The other data set was 
obtained using the COSY- 11 detection system at the COSY-Jiilich accelerator for 
Q = 3.2 MeV, 6.1 MeV and 9.2 MeV^I^. These data show enhancements of more 
than an order of magnitude in comparison to the theoretical plane wave total cross 
sections, suggesting thereby a strong FSI effect. 



2. Interpreting the near threshold pd — > pdrj data 

At low energies, ignoring the p — d FSI, an averaged squared rj ^ d production 
amplitude, (which is the ratio of the cross section for the production of the 

rj — d system to arbitrarily normalized phase space) is tradit ion ally extracted from 
data in terms of the scattering length v4,, ^ by parametrizing it as, 



where, Jb is considered a (basically unknown) scale factor and in principle depends 
upon the particular reaction studied. At the low energies considered here, this factor 
is practically energy independent, k is the relative momentum of the rj — d system. 
Within this parametrization, the analysis of the WASA data (for Q > 14 MeV), 
using different values of Aj^d taken from Ref. 20 suggests a preference for Ar,d = 
1.64 -I- i2.99 fm (as mentioned in Ref. 17). This value of the r] — d scattering length 
is obtained in Ref. 20 using an 77 — scattering length of, a^^r = 0.55 -I- i0.30 fm. 

Conclusions drawn about A^/d and subsequent inferences made about the pos- 
sible existence of quasi-bound ij — d states, from such studies, may suffer from 
limitations and uncertainties because of the reasons listed in the following two sub- 
sections. 



\m\' 
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Fig. 1. Comparison of |i^(fc)p obtained using: (a) the factorized enhancement factor prescription 
of \F{k)\'^ and (b) with the |F*'^(fc) p theoretically obtained from the full calculations generating the 
ry—d scattering wave function solving the Lippmann-Schwinger equation and using T,j constructed 
from A^^. This calculation is done for (J ~ 6 MeV. 

2.1. The scattering length description of \F{k)\^ could be 
misleading 

The scattering length enhancement factor description involves the factorization of 
the pd pdr) reaction amplitude into the 77 production amplitude and an FSI 
factor. This may not be correct. Fig. 1 exhibits this. We plot here using, (i) 

the factorized enhancement factor prescription (Eq. ^) and (ii) using the definition, 

|^th(,)|. ^ [da / dMU----^- ^ 
phase space 

where, (dcr / dAf,,^)"'™''""''''', is evaluated in a theoretical model involving the FSI 

91 

of the p, d and 77 . Here the rj ~ d FSI is included by numerically solving the 
Lippmann-Schwinger (LS) equation (as described in the next section). The t-matrix 
for rj-d elastic scattering, T,id, entering the LS equation is constructed from Ajjd 
given in Ref. 20. For the definition of (da / dAf^d)*''"™''*"^''' and the phase space 
details, see0. 

We plot the squared amplitude obtained from Eqs. (1) and (2) in Fig. 1(a) and 
1(b) respectively, as a function of the excess energy, Qnd, m the rj — d centre of mass 
system. In order to highlight the slopes of these results for different at , we choose 
values for the normalizing constant, fs, in case of the amplitude obtained from Eq. 
(1), such, that the \F{k)\'^ has the same value for the maximum value of Q,,^ (~ 3.1 
MeV) shown in Fig. 1(a). To do the same for the results obtained from Eq. (2), we 
multiply the results by a normalizing constant, TV. We observe that, as we approach 

'^For the definition of the cross section and the phase space, see Eq. (1) and Eq. (25) of Ref. 21 
and Eq. (4) of Ref. 21 or Eq. ((Sjl of this article for the t-matrix which has been used to calculate 
the cross sections. 
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the threshold, the curves corresponding to the fuU calculations for different values 
of Arid (shown in Fig. 1(b)) open up much faster than the approximate calculations 
(Fig. 1(a)). This means that the value of A^^d extracted by fitting the data by using 
Eq. (1) could be larger than that obtained by using the full calculation. 
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Fig. 2. Diflferent contributions in the inverse of the s— wave Jost function squared, [Jo{kpd)\~'^ ■ 
The eflfect of including Coulomb repulsion is also shown. 



2.2. p — d FSI is stronger at low energies 

The pd pdrj reaction, in the final state, in addition to the j] ~ d interaction, 
involves the p — d interaction, which at low energies becomes important. There are 
two possible s— wave p—d final states, the spin-doublet (^81/2) and the spin-quartet 
{'^33/2) states. ^Si/2 bears connection to the p—d distribution in the ^He bound 
state. Furthermore, since both the proton and deuteron are positively charged, at 
low energies there exists a sizable Coulomb repulsion between them too. The p — d 
FSI in t he p resent work is included through the multiplication of an inverse Jost 
function This function is known to give the modification of the non- interacting 
wave function of two particles at short distances due to the interaction between them 
(see next section for the Jost function expressions or Ref. 21 for more details). In 
Fig. 2 we show the inverse of the squared Jost function, [J{kpd)]~^ for the s— wave 
p — d interaction. In order to show the effect of the Coulomb repulsion, we have 
plotted the Jost function for the doublet and the quadruplet states with and without 
inclusion of the Coulomb part separately. The contribution of the Jost function 
without Coulomb part can be obtained by keeping = 1 in Eq. (fT3|) . We also 
show the doublet Jost function multiphed by a factor 1 -I- {\Eb\/E) to reflect the 
fact that the ^51/2 p — d state has a bound state in ('^He) with a separation energy 
\Eb\ of 5.48 MeV. We see that at small values of pd relative momentum the strong 
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interaction in the p ~ d system can enhance the wave function and hence the cross 
section, by a large amount. The Coulomb interaction, being repulsive reduces this 
effect. However, the overall effect is still an enhancement (up to a factor of 3 - 6). 

The above discussion thus makes it clear that to learn reliably about the rj — d 
interaction near threshold the data on the pd — > pdij reaction in this energy range 
should be analyzed, as far as possible, with the full theoretical framework available 
and should include all the interactions. In what follows, we have made an attempt 
to explain the new (preliminary) data available from Ref. 19 by using a formalism 

which includes the interaction between all the pairs in the final state. 

In the following section we describe the formalism only briefly and move on to 
discussing results and conclus ions. For more details of the formalism we refer the 
reader to our previous articles ^^lESES] 



3. r) Production and r) — d f-Matrices 

The cross sections for the pd pd-q reaction close to threshold are evaluated using 
the two step model for rj production, where the proton interacts with a proton (or 
neutron) in the deuteron in the first step. This produces a deuteron and a tt'^ (or 7r°). 
The pion in the second step interacts with another nucleon in the original deuteron 
and produces an rj. The two steps thus are pp{n) — > dir^ (tt*^) and vriV — > r/N. 
The ^-matrix for this model in the plane wave approximation is written as 

Tpd-^pd,j I ) = 2 * X! J (?"^ I ^) ( I Tpp(n)^d-j,+ (ttO) I ) 

I), (3) 



m s 



— mi + i e 



where fc^ is the four momentum of the intermediate pion. The summation is over 
internal spin projections and the matrix element {pn\d) represents the deut eron 
wave function in momentum space for which we use the Paris parametrization ^231 

The above kind of model for the production amplitude which involves a two step 
process (and hence a three body mechanism for production) has often been used 
for meson producing reactions at intermediate and low energies. The imp ort ance 
of such a model for the pd ^ ^He rj reaction was first shown by Laget and 
then tested with more recent data in Ref. 22 and 23. It has also been used for other 
meson producing reactions in Ref. 26 and 27. Though one worries about unitarity 
bounds on cross sections involving two step processes at high energies (77? scattering 
for example Refs. 28-30), such worries are not relevant for the present work. The 
model calculations done here are very close to threshold and in addition, the cross 
section for rj production saturates at high energies. 

The pion production in the first step being at high energy, is described by an 
on- shell p N dn t— matrix, which is parametrized using the data on this reaction 

The ( I Tpp(^n)~tdii+ {tt'>) I ) in Eq- has been related to the amplitudes ( | /| ) 
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given in Ref. 31 as 



{\Tpp{n)^d7,+ {tt») I) = \l 2 ^^^'^'^ X (I/I) (4) 

y mf-,mdPfPi 

and we relate it to the cross section using 
da mlmq \p}\ 1 

spin projections 

where pi (pf) is the initial (final) momentum in the pp ~f di: centre of mass system 
and ^JSpp-^dTT is the total energy of the same. 

The second step tt N ^ rj N being near threshold, is described by an off-shell 
matrix which is taken from a coupled channel calculation-! which reproduces the 
TT N r] N data well. We refer the reader to Ref. 32 for some recent interesting 
investigations of the vr N data. 

To take the FSI into account we include the interactions between the p — d and 
r]~d pairs. The rjp interaction is implicit in our description of the production vertex, 
which uses the matrix for tt N — > 77 A'^. Though, in principle, a finer treatment 
of the three body problem (involving the 3-body p-d-rj scattering) would be more 
appropriate, we expect the present approach to work reasonably well for the problem 
under consideration. Some remarks in connection with this issue are in order here. 

The choice of any approximation is guided by the characteristics of the problem 
at hand. In our case, because of the large momentum transfer (^ 770 MeV/c), the 
reaction takes place over a small volume of space. In such situations (first intro- 
duced by Watson and M igdal and subsequently discussed in most of the scattering 
theory text-books ^^^^^H)^ an essential feature of the cross sections is that the energy 
dependence of the cross section due to FSI, factors out from that of the primary 
production amplitude. The FSI factor, which has the correct high energy behaviour 
is usually taken to be the inverse of the Jost function. In the present paper we have 
used this prescription for the p-d interaction. Having done that, we could include 
the eta-d interaction fully, hence we did not resort to a similar factorization for it. 
We included the eta-d interaction effect fully by keeping it inside the production 
amplitude and using an exact scattering wave function. The above factorization 
method hasbeen used in literatu re to describe the low momentum behaviour of 
the meson ESIMI and pion spectra [ 3'^l'^^l measured in proton-proton collisions, and 
study the total cross sections in the — > p A'+ A reaction 1221, 

The transition matrix for the reaction, pd — s- pdrj, which includes the interac- 
tion between the rj meson and the deuteron is written as, 

T= {'^]^^^, k'j,]m'p,m'^\Tpd^pdri\kp,kd]mp,md) (6) 

where kp and kpi are the proton momenta in the initial and final states, respectively. 
krid is the 77 momentum in the i] — d centre-of-mass. , , mp' and m^' are the 
spin projections for the proton and the deuteron in the initial and final states. The 
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final state 77-d wave function, * consists of a plane wave and a scattered wave 
and satisfies the Lippmann Schwinger equation for ry — d elastic scattering: 

Here, T^d is the f-matrix for rj-d elastic scattering. Replacing the above equation in 
the transition matrix for the pd — > pdrj reaction (Eq. ([6])), leads to 

T — {kjjd : kpi ; rrLpi , m^' \ Tp d^pdril kp , kd] , irid ) 
dq {kr,d\ md- \ T,^d\q; rn2' ) 



E 



(27r)3 Eikr,d) - E{q) + le 



x{q, kp, ; m2' , mp/ |Tpd^prf^ \ kp ,kd; nip , rud) , (8) 

where the first term represents the ry production without any interaction with the 
deuteron and the second one takes care of the rescattering of the 77-d pair to all 
orders. It is indeed the half-off-shell t- matrix, Trjd, which is reponsible for converting 
the off-shell eta mesons in the intermediate state to the on-shell ones in the final 
state. 

Since the low energy r/N interaction is dominated by the Sn resonance 
N*(1535), we perform a partial wave expansion of the rj-d i-matrix, Trid{krid-iQ) 
(= [krid] nid' \ Trid\q', ni2')) and retain only s-waves. Trj d{kri d, <t) is half off-shell 
and we write it as a product of an on-shell t— matrix and an off-shell (which is 
actually half off-shell here) form factor. Thus, we write Trjd{krid,(t) using s-waves 
only as, 

T^d(k^d,q) = T°^d{k^d) x .g(fc,<j, g), (9) 
where the T°^(krjd) is written in terms of the 77c? scattering length as 

r,rj--^x(-L_,fc,,A (10) 

We choose to define the off-shell form factor in terms of the deuteron form factor 
as in Ref. 40, 



gikr,d,,q)^ J dfjo{rk^d/2)\Mr)r.Mrq/2), (11) 

where an average over the directions of the momenta has been performed in order 
to extract the s-wave contribut ion. The deuteron wave function, 0rf(r), is written 
using the Paris parametrization ^241 

The r/d scattering length, used in Eq. (jlOp. is taken from Ref. 20, where it 
has been deduced by solving three-body Faddeev type equations. This calculation 
includes the scattering of the 77 on the break-up continuum states as well as on the 
deuteron in the intermediate state. We use the A.,, ^ in this work calculated for a 
wide range of rj — N scattering lengths since the rj — N interaction is still not very 
well known. 
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For the d interaction, as mentioned above, we use the Jost function, which is 
justified because the rj production involves high momentum transfer. Near threshold 
this value is around 700 MeV/c. This means that the reaction occurs in a small 
interaction volume. The Jost function is given as in Ref. 21, 

\Eb\. 



D 



(12) 



The expressions for spin quadruplet (Q) and doublet (D) [Jo {kpd)] ^ are given by 



[Jo {kpd)\x — 



A 



X 



1 + cot^ 5x 



(13) 



where X is either Q or D, = 1 and Aq = 2. The factor above originates 
from the Coulomb interaction. Details regarding the parameters used in Eq. 
can be found in Ref. 21. 

The final amplitude is calculated by multiplying Eq. ([12]) to the squared i-matrix 
obtained after including the rj — d interaction, i.e., Eq. ([8]). 



4. pd and r/d FSI effects on the total cross sections 

We present the total cross sections calculated, using our full amplitude described 
above, along with the data for excess energies up to 20 MeV in Fig. 3. The COSY- 
11 data shown in Fig. 3, though preliminary, seems to be in reasonable agreement 
with the previous data near threshold. In this figure we show separately the relative 
contribution of the different pieces in the FSI. The rj — d FSI has been calculated 
using an ?7 — c? scattering length corresponding to o^at = 0.98 + iO.37 fm. We find 
that (i) both, the rj — d and p — d FSI affect the results equally strongly and in 
the whole region of the excess energies, (ii) Comparison of results with and without 
Coulomb repulsion show that the Coulomb interaction reduces the cross section by 
a large amount. In fact, it is because of this interaction that the calculated cross 
sections come in accord with the experiments. This is one of the main result of our 
work. 

We have also studied the sensitivity of the results to the rj — N scattering length 
which is still not so well known. The results shown in Fig. 4(a) have been obtained 
using the Paris parametrization for the deuteron form factor and four different 
values of the rj — d scattering lengths, A,, rf. These A,, ^ are calculated for different 
input values of ar/N, ranging from a weak to a strong rj — N interaction. For the 
results presented here, these a^Ar values are 0.25 + i0.16, 0.44 + i.30, 0.579 + ?0.399, 
and 0.98 + iO.37 fm. The different A^jd results in the figure are identified by their 
corresponding a^jN values. We see that (i) all the calculated cross sections show a 
strong FSI effect, and (n) the results for a^Ar = 0.44 + i.30 and 0.579 + zO.399 fm 
seem to pass through the maximum number of data points. 

For the sake of completeness, in Fig. 4(b), we show the results obtained by 
using in Eq. (fTT|) , prescriptions of the deuteron wave function other than the Paris 
form, which are also sometimes used in literature. We have thus calculated the total 
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Fig. 3. Relative contribution from different pieces in the FSI. The results are for a^ pf = 0.98 - 
jO.37 fm. 
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Fig. 4. A comparison of data | 17 | 18 | 19 | ^j^]-, ^^j^g total cross sections calculated including all the 
final state interactions, (a) using the Paris parametrized wave function for the deutoron in r? — d 
FSI and varying the values of a^jv and (b) for a^ n = 0.44 + i.30 but using different deuteron 
form factors in the FSI. 



cross sections using the wave functions obtained from the Paris potential, Hulthen 
potential, and the Gaussian form. The scattering length ariN for the calculation in 
Fig. 4(b) has been kept fixed at 0.44 + «.30 fm. The maximum difference in the total 
cross sections, obtained due to the use of different wave functions, is found to be a 
factor of the order of 1.6 at higher excess energies. 
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5. Summary 

We have presented a theoretical analysis of the data on the pd pdij reaction near 
threshold using a model which was reasonably successful in reproducing the data at 

higher energies. Wc start by comparing the model calculations with a form of the 
scattering amplitude often used for data analysis and investigate the limitations 
of this form. The cross sections for the pd — > pdr] reaction at low energies (up 
to ~ 20 MeV above threshold) are then evaluated within a two step model for rj 
production to investigate the importance of including the FSI between all particles 
at low energies. The conclusions of this study can be summarized as follows: 

(1) Quantitative reproduction of the experimentally observed enhancement in the 
cross section near threshold, requires the inclusion of the interaction between 
all three particles, namely, the proton, deuteron and the eta meson in the final 
state. Both the strong and Coulomb p — d FSI are found to be important at low 
energies. 

(2) The best agreement with experimental data seems to come from the //-nucleon 

scattering length values of a,, tv = 0.44 + i.30 fm and 0.579 + i0.399 fm. Calcu- 
lations based on the exact Alt-Grassberger-Sandhas equations in Ref. 20 lead 
to T] — d scattering lengths of (1.15 + il.89) fm and (0.34 + i3.31) fm for the 
two above Ujjn's respectively. 

(3) In the energy range discussed here, the sc;attering length description of the 
T] — d ^-matrix, T^d, which enters the Lippmann-Scliwinger equation for the 
T] — d elastic scattering wave function, seems to be sufficient. 

(4) The extracted value of the ?7-deuteron scattering length, A,, d, from experimental 
|F(A;)p using the simplified relation in Eq. (1) can be misleading. 
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